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We extend our previous work on soliton ratchet devices [L. Morales-Molina et al, Eur. Phys. 
J. B 37, 79 (2004)] to consider the joint effect of two ac forces including non-harmonic drivings, 
as proposed for particle ratchets by Savele'v et al. [Europhys. Lett. 67, 179 (2004); Phys. Rev. 
E 70 066109 (2004)]. Current reversals due to the interplay between the phases, frequencies and 
amplitudes of the harmonics are obtained. An analysis of the effect of the damping coefficient 
on the dynamics is presented. We show that solitons give rise to non-trivial differences in the 
phenomenology reported for particle systems that arise from their extended character. A comparison 
with soliton ratchets in homogeneous systems with biharmonic forces is also presented. This ratchet 
device may be an ideal candidate for Josephson junction ratchets with intrinsic large damping. 
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The understanding of ratchet mechanisms is a very ac- 
tive field, of wide interest by its potential application 
in the design of devices with new transport properties. 
The key feature of ratchet devices is their ability to rec- 
tify the motion of particles subjected to an external ac 
force with zero time average. Originally proposed as a 
toy model for molecular motors, in the last decade many 
proposals have been put forward for devices that use this 
ratchet phenomenon in different applications Ratch- 
ets working with extended particles (solitons) were sub- 
sequently studied as a generalization of particle ratchets 
0- Recently, a ratchet device driven by two ac forces 
was introduced by Savele'v and coworkers Q, in which 
the combination of the two drivings produced a variety 
of interesting phenomena and allowed a finely tunable 
control. In view of the rich behavior demonstrated by 
this system, a very natural issue is its extension to soli- 
ton ratchets, that have important applications which can 
benefit from this proposal. 

Soliton ratchets under the presence of two ac forces 
have been extensively studied 0, 0, El m homogeneous 
systems, i.e., without an underlying ratchet potential. In 
this case, it has been shown that the ratchet mechanism 
works only for asymmetric biharmonic forces, where the 
appearance of a directed translational motion is a result 
of the effective coupling between the internal mode (os- 
cillations of the soliton width) and the external driving 
force. In this paper, we focus on a soliton ratchet device 
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recently studied by us 0, |j| , based on a nonlinear Klein- 
Gordon model with a lattice of delta-like inhomogeneities 
that induce a ratchet potential for the solitons. As we 
will see, one of the advantages of our model as compared 
to the homogeneous system is that it works irrespectively 
of the symmetry of the ac driving: Directional transport 
takes place for ac forces with commensurate frequencies. 
Moreover, in the homogeneous system the motion drasti- 
cally decreases for increasing damping, due to the slow- 
ing down of the soliton width oscillations , while in the 
present system the ratchet phenomenon is present up to 
rather large values for the damping. 

Our model, first introduced in |8|, is given by 

4>tt + P4t ~ <f>x* + bLd(0)[1 + V(x)] = f(t), (1) 

where fit) = A\ sin(wit + Oi) + A2 sin^i + #2), being 
Ai, A 2 the respective amplitudes of forces, ui, UJ2 the 
frequencies and 6*i and 62 the phases of the harmonics 
(see e.g. Q). Here we focus on the sine-Gordon model, 
although the same scheme can also be applied to the 
general nonlinear Klein- Gordon model 0. For V(x), we 
choose a spatially periodic potential, where the unit cell is 
given by an asymmetric array of delta functions (inhomo- 
geneities) in order to produce a ratchet-like phenomenon. 
Specifically, the unit cell, of length L, is defined by three 
inhomogeneities with the same intensity, the first one lo- 
cated at the beginning of the cell, the second one at a 
distance a from the first one, and the third one at a dis- 
tance b from the second one, i.e., 

V(x) = e [S(x — xi — nL) + S(x — X2 — nL) 

n 

+S(x -x 3 - nL)} , (2) 
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FIG. 1: Mean velocity vs damping coefficient for different 
frequencies, squares: u — 0.015; filled circles: uj — 0.05. The 
parameters used are e = 0.5, A = 0.2, x\ = 0.5, x 2 = 1, 
x- A = 2.3 and L = 4. 



where L = a + 6 + c, a = 2:2 — ari, b — X3 ~ X2 and 
c = L + xi — X3, with xi < a: 2 < X3. The parameters 
(a, b, c) are chosen to be comparable to the static soliton 
width Iq in absence of inhomogeneities. In addition, these 
should fulfill the conditions a,b < c with a ^ b. 

As a preliminary step, we begin our study by analyz- 
ing how the system works in different regimes of damp- 
ing. This is very relevant to applications such as Joseph- 
son junctions. While the standard Josephson junctions 
work usually at small damping, junctions with intrin- 
sically high damping such as superconductor insulator 
normal-conductor insulator superconductor LJJ or high- 
T c LJJ technology can also be fabricated 0. For this 
particular aspect, we look at the case of a single har- 
monic component in Eq. In this case, according to 
0, |i| we have a directed motion whose direction is de- 
termined by the position of the inhomogeneities, and the 
dynamics reduces to a system very much like a rocking 
ratchet for a single particle [Tol ] . Our results are shown in 
Fig. ^ where we have chosen a lattice of inhomogeneities 
whose configuration (a < b) yields a negative current @ . 
Hereafter (X) means the time average of the velocity. 
The figure demonstrates the increment of the efficiency 
upon decreasing the damping /3; however, nonzero values 
for the velocity are found for rather large damping val- 
ues, and for the lower frequency there is still motion even 
for 0>1. 

We note also that for these frequencies with lower 
damping values, the dynamics results depend on the ini- 
tial conditions, which may lead to a chaotic dynamics. 

Let us now move to the situation with two simultane- 
ously acting ac forces. For the time being, we work with 
sinusoidal forces as in Eq. (JTJ , and we will consider later 
the case of rectangular pulses as in |j. With two har- 
monics present (doubly rocked ratchet) we have a system 
in which the symmetry can be broken both spatially (re- 
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FIG. 2: Mean velocity vs force amplitude A = A\ 
Simulations of Eq. the parameters used are = 1, o»i = 
0.025, L0 2 = 0.3, e = 0.5, 61 = 2 = 0, xi = 0.5, x 2 = 1, 
xs = 2.3 and L — 4. The solid line is a guide for the eye. 

flection symmetry V{x) = V{—x)) and temporally (time 
shift symmetry f(t) = -f(t + T/2) with T = 2tt/uj{). 
As we will now see, it is possible to obtain current rever- 
sals irrespective of the symmetry of the biharmonic force. 
The reason is that the underlying mechanism of trans- 
port studied for homogeneous systems subjected to bihar- 
monic forces, that required temporal symmetry breaking 
[6j , is not responsible for the directional transport in our 
present model, as can be inferred from the fact that in the 
inhomogeneous system we have directed transport under 
large damping. We examine the dynamics of Eq. Q in 
the partially adiabatic regime u>i/u>2 <C 1 with u>\ <C 1 
and u>2 < 1, setting A\ = A%. In this case, multiple cur- 
rent reversals are possible in the particle ratchet system 
Q. Figure El confirms clearly the existence of several cur- 
rent reversals in our soliton device. The results show that 
we can reverse the direction of motion by changing the 
force amplitudes, which, as in the particle case, opens the 
possibility to control the rectification properties in great 
detail. 

In the previous analysis we have taken /3 = 1, i.e., a 
rather large damping where the inertial effects are small. 
However, from the above discussion for one harmonic, 
we expect changes in the behavior also in the case of a 
biharmonic force when the damping is reduced. Figure 
13 exhibits indeed a different picture for the dynamics as 
compared with Fig. [3 The main differences are related 
to the shift of the windows of motion towards lower force 
amplitude values as well as an increase of the absolute 
value of the mean velocity for some force amplitudes. 
Interestingly, we note from this picture a non-typical be- 
havior for the average velocity close to the region where 
the first current reversal takes place: We observe that, 
while for some values of the force amplitude the windows 
of motion are suppressed, for other cases, the absolute 
value of the average velocity is enhanced and even re- 
versed. 
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FIG. 3: Mean velocity vs force amplitude. /3 = 0.4, A\ = 
A2 — A simulations of Eq. The other parameters are the 
same as in Fig[5] The solid line is a guide for the eye. 
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FIG. 4: Mean velocity vs 6. Here we have taken 9i = and 
82 = 6 in Eq. The other parameters are A\ — A2 = 0.16, 
P = 0.4, u)i = 0.025 and lu 2 = 0.3. 



This issue can be further examined in Fig. 0] where the 
dependence of the average velocity on the relative phase 
of the drivings is plotted. Here we see that by chang- 
ing the phases we reverse the direction of motion. This 
behavior is not only restricted to this singular region in 
FigED m i'ACt, it can be observed close to the regions 
where the currents are reversed in FigEI We stress that 
we have taken 8\ = and 82 — 8 without loss of general- 
ity, since although the dynamics obviously changes with 
both phases 8\ and 82, one can always map the choice 
for 81 and 82 into the previous representation through 
the transformations 8' 2 — 82 — (^2/^1)^1 with 8[ — 0. A 
feature of the behavior for this doubly rocked ratchet sys- 
tem, which differs from the homogenous system case, is 
the quantization of the velocity dependence on the phase 
(see Fig. 01 and cf. Fig. 2 in 6]). This is yet another 
evidence that the ratchet mechanisms at work in the ho- 
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FIG. 5: (Color online) Mean velocity for different frequency 
ratios 0^2/^1 in the fully adiabatic regime with rectangu- 
lar wave signals f(t) = Ai(t) + A 2 (t). The parameters are 
wi = 0.005, Ai = 0.25, A 2 = 0.05, /3 = 1 and 9i = 2 = 0. 
Inset: Rectangular wave signals: Ai(t) = Ai sgn[sia(u)it+9i)] 
(solid line); Ai(t) = A2 sgn[sin(tJ2t + (dashed line); 
A 2 = A-l/2 = 0.1, w 2 = 2u>i, 9i=0 and 2 = it/2. 



mogeneous and inhomogeneous cases are not the same. 
Nevertheless, despite this quantized nature, it is feasible 
to control the direction of motion by tuning the phases 
of the biharmonic force as in the homogeneous case. 

In view of the similarities with the single particle 
ratchet system in it is important to assess the de- 
gree of similarity between the two systems. Savele'v and 
coworkers focus mainly on the case of rectangular wave 
signals as the asymmetry and nonlinearity-induced mix- 
ing are separable ■ For such rectangular wave signals in 
the fully adiabatic regime (0^1,072 "C 1) with a sawtooth 
ratchet potential, Savele'v et al. report changes in the dy- 
namics only for a relation between the frequencies given 
by uj%/wi = (2m — l)/(2n — 1). As Fig. ^ shows, in our 
case the behavior turns out to be much more complicated. 
The average velocity fluctuates around the value indi- 
cated by a horizontal line, but it does exhibit very many 
peaks. Most importantly, unlike the single particle pic- 
ture shown in Q , here the peaks appear not only for frac- 
tional harmonics (^2/^1 = (2m— l)/(2n— 1), but also for 
harmonics that fulfill the relations = (2m — l)/2n 

and = 2m/(2n — 1). 

The reason for the difference between the particle and 
the soliton doubly rocked ratchets can be traced back to 
the extended character of our solitons. It is well known 
that the soliton dynamics is affected by the deformation 
of the soliton width that accompanies its motion along 
an array of inhomogeneities |9| . This in turn modifies the 
effective potential arising in the description of the soliton 
as point particle, as well as in the variations of the cor- 
responding effective force. Accordingly, there is a large 
degree of feedback between the soliton width and the soli- 
ton motion. To obtain some insight on these issues, it is 
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FIG. 6: (Color online) Collective Coordinates: Mean veloc- 
ity for different frequency ratios 012/ loi in the fully adiabatic 
regime. u>\ = 0.005. The other parameters are A\ = 0.25, 
A 2 = 0.05, f3 = 1 and 0i = 62 = 0. 



necessary to study the evolution of the degrees of freedom 
that are involved. In order to do so, we resort to the use 
of the collective coordinate approach ^lj which involves 
the soliton width as an additional degree of freedom. In 
doing so, we find that the soliton can be described by two 
collective variables X and I, whose expressions are given 
by Eqs. (4)-(5) in 0, with the difference that now f(t) 
contains the two rectangular wave signals whose expres- 
sion appears in the caption of Fig. El This is in contrast 
to the single particle doubly rocked ratchet of Savele'v 
et al, because of the appearance of the width degree of 
freedom. This collective coordinate approach explains 
in full detail the results of the simulations, as shown in 
Fig. for the fully adiabatic regime. Again, as in the 
simulations of Eq. (JIJ we note not only the existence of 
peaks for the frequency ratios 102/^1 = (2m— l)/(2n— 1), 
but also for some frequencies ratio that fulfill the relation 
LJ2/0J1 = (2m — l)/2n and W2/0J1 = 2m/ (2n — 1), absent 
in the single particle picture y|. This result is a clear 



demonstration of the role of the soliton width in the dy- 
namics. 

To conclude, we have generalized the results obtained 
for doubly rocked particle ratchets to extended sys- 
tems, finding new phenomena that arise from the intrin- 
sic width of the solitons. It was shown that the direction 
of motion can be modified by changing the relation be- 
tween the frequencies ratio, the phases of the harmonic 
forces as well as their amplitudes. However, in the frame 
of the fully adiabatic regime of our soliton ratchet we find 
many more peaks than in the particle ratchet system of 
Savele'v et al, as we see special velocities for several types 
of frequency ratios. We have been able to show that this 
is thoroughly accounted for by the interplay of the soli- 
ton width and motion. We have also compared to doubl 
rocked soliton ratchets in homogeneous systems 0, |T 
and verified that, although the soliton width is involved 
in both cases, the mechanism for the appearance of the 
ratchet effect is different. Aside from the fact that homo- 
geneous soliton ratchets arise only for asymmetric bihar- 
monic drivings, further important differences include the 
damping dependence of the velocity and the quantization 
of the dependence of the velocity on the relative phase. 
We emphasize that our ratchet system can be straightfor- 
wardly implemented in a Josephson junction device 0. 
In that case, the very many possibilities for the motion 
we have reported here would allow for a highly control- 
lable device that can be tailored to fit different specific 
applications. The property that the ratchet phenomenon 
is observed even in the presence of large damping makes 
this system preferable to a homogeneous one driven by 
a biharmonic force, and makes it much more suitable for 
real life applications. Finally, we note that an analysis 
similar to the present one can be extended to incommen- 
surate ac forces with irrational values for the relation 
LJ2/LO1. However, our preliminary results show that this 
choice gives rise to new phenomena that require careful 
attention, and therefore it will be the object of a future 
investigation. 
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